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ABSTRACT 

We show that the inter-cloud Larson scaling relation between mean volume density 
and size p oc R , which in turn implies that mass M cx R 2 , or that the column 
density N is constant, is an artifact of the observational methods used. Specifically, 
setting the column density threshold near or above the peak of the column density 
probability distribution function iV-pdf (N ~ 10 21 cm~ 2 ) produces the Larson scaling 
as long as the 7V-pdf decreases rapidly at higher column densities. We argue that the 
physical reasons behind local clouds to have this behavior are that (1) this peak column 
density is near the value required to shield CO from photodissociation in the solar 
neighborhood, and (2) gas at higher column densities is rare because it is susceptible 
to gravitational collapse into much smaller structures in specific small regions of the 
cloud. Similarly, we also use previous results to show that if instead a threshold is 
set for the volume density, the density will appear to be constant, implying thus that 
M oc R 3 . Thus, the Larson scaling relation does not provide much information on 
the structure of molecular clouds, and does not imply either that clouds are in Virial 
equilibrium, or have a universal structure. We also show that the slope of the M — R 
curve for a single cloud, which transitions from ncar-to-flat values for large radii to 
a = 2 as a limiting case for small radii, depends on the properties of the iV-pdf. 
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1 INTRODUCTION 



More than 30 years ago, iLarsonl (| 198 lT ) published his scaling 
relations for molecular clouds (MCs): the column density- 
size relation (or "Larson's third relation") 



pociT 1 

and the velocity dispersion-size relation 
dv oc R 12 . 



(1) 



(2) 



272 — 2 will be consistent with V irial equilibrium (see e.g., 
IVazauez-Semadeni fe Gazollll995l h 

Larson's relations have been used in many pa- 
pers t o describe the inte rnal structure of clouds 
(e.g., iGoldbaum et aL 2011, to cite just one of the 
more recent examples). However, their validity has 
been called into quest i on for observational reasons 
jKegell Il989t IScald fl990l: IVazquez-Semadeni et~afl Il997l; 



The exponents reported in that work were 71 = —1.1 
and 72 = 0.39, respectively. However, the more widely 
accepted values are 71 = — 1 and 72 = 0.5; and these 
have been usually thought to be observational evidence 
that clouds are in Virial equilibrium, although it is im- 
portant to recall that any pair (71,72) satisfying 71 = 



* e-mail: j .ballesteros@crya.unam.mx 



Ballesteros-Paredes fc Mac Low! 2002: iBallesteros-ParedesI 
20061 : iBallesteros-Paredes et alj|2011a ). Specifically, clouds 
should have a minimum column density to be detected; and 
if the column density where too large, clouds will become 
optically thick (in 12 CO), and then such dense regions 
would not be easily detected. The net result of these effects 
would be that observed CO clouds should exhibit a small 
dynamic range in column density, independent of their 
intrinsic structure. 

In a recent study, lLombardi et al.1 (|2010l ) used near- 
infrared extinction measurements to increase the dynamic 
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range of inferred mass column densities. They studied a 
sample of local MCs with substantially different physical 
properties, from giant MCs hosting massive star formation, 
and masses of the order of 2xl0 5 Mefe.g., Orion), down 
to clouds with less star formation activity, and masses of 
the order of several thousand Mq (the Pipe Nebula). Their 
main results are: 

• For an ensemble of MCs observed at a given extinction 
threshold Ao, the mass varies with size as M oc R 2 , imply- 
ing that the column density E = N '/ ' \xrriH is constant and 
depends on that threshold. 

• The numerical value of the column density changes with 
the value of the threshold Ao. 

• When analyzing a single MC at different thresholds Ao, 
the mass-size relation varies as a power-law, with a some- 
what shallower power law index {M oc r l e ). 



lLombardi et~ai1 l|2010l ) argue that their results can be 
understood as a consequence of the fact that clouds ex- 
hibit a column density probability distribution function (N- 
pdf) w ith a lognormal function al form, alth ough in a recent 
paper , iBeaumont et al.l l|2012l ) revisit the lLombardi et al.l 
(2010) analysis, and suggest (without demonstration) that 
other forms than a lognormal iV-pdf can produce the same 
result. They also argue that the only requirement to satisfy 
Larson's third relation is that the mean iV-pdf is uncorre- 
cted with region-to-region dispersion in area. 

In the present work we discuss why MCs with very dif- 
ferent properties appear to have very similar column densi- 
ties at a given threshold. In fj2] we dr aw S — Aq diagrams to 
show that the main features found bv lLombardi et all (|20ld ) 
are found for very different functional forms, and that the 
scatter in column density at any threshold is smaller than a 
factor of 3. In i|4]we show, furthermore, that when thresh- 
olding the volume density, what appears to be constant is 
the volume density, suggesting that the relation is a conse- 
quence of the procedure adopted. We argue that the only 
special feature of the structure of local MCs is that they 
peak at a particular value of the mean iV-pdf, in agreement 
with IBeaumont et al.1 ||2012T ). In §5\ we analyze the case of 
the mass-size relation for a single cloud, showing that rather 
than a power-law, the relation is a curve that has slope ap- 
proaching to 2 at small radii, and flattens for larger radius. 
As a consequence, there is some interval in which the slope 
is close to 1.6-1.7, as reported in observations. The proper- 
ties of this relation are determined directly by the shape of 
the iV-pdf. In iJ5]we further identify the physical reasons for 
having a typical iV-pdf: First, there is a minimum column 
density for CO shielding; and second, above this threshold 
the amount of mass decreases rapidly with increasing col- 
umn density, which is probably the result of gravitational 
contraction and collapse. 



iKainulainen et all |2009) have a peak at Ah ~ 0.1, equiva- 
lent to Ay ~ 1 or N ~ 10 21 cm -2 . Now, given a probability 
distribution function of the extinction p(A)dA, the area and 
mass above some extinction threshold Aq are defined as: 



(3) 



S(Ao) = Stot / P(A) dA 

JAo 

and 

r °° 

M(A ) = Stot^mn A p(A) dA 



(4) 



where \i is the mean molecular weight within the MC, itih is 
the mass of the hydrogen atom, and ft is the ratio between 
the total h ydrogen column density and the dust extinction 
coefficient. lLombardi et all l|2010l . see their Fig 4) showed 
that if clouds have a lognormal iV-pdf, the surface density 
as a function of Ao, which is given by the ratio of eq. Q 
divided by eq. <(3j , is nearly constant for low thresholds Ao, 
and varies approximately linearly for large thresholds, i.e., 



E = 



constant 
oc A 



if Ao « A 

max • 

if A » A max , 



(5) 



where A max is the extinction value at which the iV-pdf peaks. 
Since this behavior reproduces qualitatively that of the 12 
MCs in the Solar neighborhood they analyzed, these authors 
suggest that the observed constancy of the column density 
for different clouds is due to their intrinsic internal structure, 
given by a lognormal iV-pdf. However, many of the MCs in 
the solar neighborhood do not exhibit a lognormal iV-pdf at 
large column densities but rather fr equently exhibit power- 
law or more complicated shapes (|Froebrich et al.l 120071 ). 
Thus, the behavior of E(Ao) described above cannot be con- 
sequence of a universal distribution function of the column 
density of the MCs. 

In order to understand the origin of Larson's third rela- 
tion for multiple clouds, we analyze the following examples. 
For clarity's sake, we will first describe the examples, leaving 
the discussion for 



3 THE INTER-CLOUD MASS-SIZE RELATION 
FOR SIMPLE COLUMN DENSITY PDFS 

3.1 A single power-law 

First, assume that the iV-pdf has a power-law shap43 : 

P {A)dA = po{A/A m ^y n dA (6) 

Using this equation, the mass above some extinction thresh- 
old is given by 



M(A ) 



Stotf-irriH /3poA max / Ao 



2-n 



A n 



(7) 



FORMALISM: COLUMN DENSITY PDFS 

AND THEIR E — Ao DIAGRAM 1 Note that plotting vs. dlogj4 will produce a power law with a 



In our analysi s we follow the formalism outlined by 
lLombardi etafl (|20ld ). We use A to denote the extinc- 
tion in the K band unless otherwise specified. We also re- 
call that Ay ~ 10 Ak , and that the iV-pdf published by 



different index, 
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were we have assumed that n > 2 in order to get a finite 
mass (indeed, by simple inspection one can verify that all 
the column density histograms of MCs published in the lit- 
erature fall faster than a power- law with an exponent o f 
m ^ 2.5, i.e., n 3.5, see e.g., iKainulainen et all [2009) . 
Similarly, the area of the cloud will be given by 



S(A ) 



A 



1 - n \ A n 

and thus, the column density will be simply given by 
^_ M(A ) , , (I - ;,. \ / 1„ 



S(A ) 



fxirmPAx, 



An 



(8) 



(9) 



We note then, that for a power-law distribution, the 
column density alway s varies linearly with Ao (see also 
iBeaumont et al.ll2012h . 



3.2 Two power laws 

Now assume that the iV-pdf has a two power-law shape: one 
with a positive exponent, for low extinctions, and the other 
with a negative exponent, for large extinctions, i.e., 



p(A) 



Po(A/A max ) 9 for A < A mm , 
p (A/A max )" n for A ^ A max , 



(10) 



with q, n > 0. Plugging this 7V-pdf into eqs. ([3]) and Q we 
obtain 



S = Ki 



1 



,1 + 5 
with 

Ki = StotPoA, 
and 

M = Ki ' ' 



Aq 



1+8' 



1 



1 -71 



2 + q 



An 



An 



2+q- 



2-n 



with 



K2 = StotfimHPpoAn 



(11) 



(12) 



(13) 



(14) 



Making K3 — K2/K1 — iituh /3A max , the surface density is 
given by 



E(A ) = K s 



l-(A /A m ^) 2+q 


(2~ 


- 9 )- 1 -(2-n)" 1 


1-(A /A max ) 1+9 


(1- 


-g)-i-(l-n)-i 



(15) 



for Ao < imax- For larger Ao, the case is reduced to a single 
power-law, and the surface density will be given by eq. ([9]). 

In Fig. \T\ we plot E vs Ao for a two power-law A^-pdf. 
In both cases, we have assumed A max = 0.1, the typical 
val ue found in a larg e samp le of Solar neighborhood MCs 
by IKainulainen et all |2009l V In the upper panel we vary n 
from ~ 3.8 to 5.8, in steps of 0.5. This is approximately 
the range of exponents that can be found in the 7V-pdfs 
presented by IKainulainen et all (|2009l . m ~ 2.8 to 4.8.) As 
expected, the change in the power-law index modifies the 
surface density at every Ao- In the lower panel, however, 
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Figure 1. S vs Ao, from cloud 7V-pdfs given by two power laws, 
as in eq. JlOjl. The upper panel presents curves with q = and 
n = 2.8, 3.3, 3.8, 4.3, and 4.8. Lower panel: curves with n =, and 
q = 1 to 5.5 in steps of 0.5. Note that in both cases, the inferred 
column density S is the same, within a factor of ~2. 



varying the low-extinction index q (from to 4.5 in steps of 
0.5) only modifies the surface density for low Ao- 

Two more features can be inferred at this point. First 
of all, for large Ao, this case is reduced to the previous one 
(a single power-law), and then, the column density varies 
linearly with the extinction threshold, i.e., E oc Ao. On the 
other hand, at low Ao (<< i m ax), the column density for a 
cloud with a given iV-pdf becomes constant: 



E(A ) ~ /j,m H /3A n 



1 + 5 

2 + q 



(16) 



We notice, moreover, that for typical values^ of n, q, the 
minimum column density given by eq. (|16[) varies less than 
a factor of two. 



3.3 A lognormal A-pdf 

The case of a lo gnormal A-pdf was analyzed by 
lLombardi et al.l (|2010h . and the mass surface density is given 
by 



2 Typically (see IKainulainen et alj feOQQ). 3.8 ^ n ^ 5.5 and 
— 1 ^ q ^ 5 (recall that q = 1 corresponds to a flat JV-pdf, 
according to footnote [TJ , where m is the slope measured in the 
N-pdf and q plays the role of n. 
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Figure 2. S vs Ao, for a cloud that exhibits a N-pdi given by 
a lognormal function. Here we assume Amax — 0.1. The different 
curves follow eq. H17H with cr = 0.29 and 0.59 ( the range of varia- 
tion in a for the lognormal functions fitted by Kainulainc n et al.l 
2009), plus two extreme cases: a = 1, and 1.5. We note that, even 
for these unrealistic cases, the inferred column density E is the 
quite the same, within a factor of ^ 3-5. 



E(A ) 



/jmji^max exp 



erf {[In (A /A max ) - a 2 }/V2a} 
1 - erf [In (A /A max )/V5o-] 



(17) 



Although the non-dimensional plot of this function has 
been already presented and discussed by these authors, for 
comparison we display it in dimensional form, assuming 
A m « = 0.1. The qualitative behavior is similar to that of 
the iV-pdf exhibiting two power-laws previously discussed. 

We further note that, since the error function can be 
expanded as 



,/ \ , exp(-:r. ) I _ x 1 _ 3 3 _ 6 . 
erl(:r) ~ 1 — — - \x — —x + -x + 



(18) 



to first order the mass column density can be written as: 



E - /im H /3A max exp ( ^— ) ln(A /A max ) x 



exp [— (ln(Ao/A 

11 



a 2 ) 2 /V2a] 



_ (ln(A /A max ) - a 2 ) CX P [- hl( A /A majf ) /2a 2 ] 

for large values of Aq. Simplifying, 
(A /A max )ln(A /A 



HturPA^ 



In (A /A n 



(19) 



This equation makes clear that, for large Ao/A n 
ln(Ao I Arasx) » c 2 and thus, 

E oc A . 



(20) 



On the other hand, for small Ao, E becomes constant with 
a value of 



E(A ) = /xm^ ,9A max exp ■ 



(21) 



3.4 Column density PDF with several peaks 

We now proceed with a slightly more complicated example: 
a three-peaked iV-pdf defined by three lognormal functions, 
as follows. 



p(A)dA = § ex P 

i=l,3 



(In A - In Ai 
2a 2 



dA, 



(22) 



and thus, the surface density will be given by 



271 



E 



Pl A t e^ /2 [l - erf ((a, - In (Aq/AO)/^)] / 



Y Pid [l " erf [In (Ao/A*)/^]] 



(23) 



In Fig[3]we show our hypothetical iV-pdf (upper panel), 
and its respective E vs Ao function (lower panel). The pa- 
rameters of the lognormal functions were chosen to enhance 
their effect on the column density. In reality, the observed 
iV-pdfs exhibit smaller variations. As can be noted from this 
figure, we again reproduce the main features observed in the 
E — Ao diagram: for Ao >> Ai, E increases linearly with Ao, 
and for for Ao << Ai, E becomes constant. Additionally, in 
this case we note that the peaks in the 7V-pdf produce small 
fluctuations with a trend that has a slope close, but smaller 
than unity. We will get back to such fluctuations below. 



3.5 Lognormal with a power-law tail PDF 

Most of the observed column density histograms appear to 
behave as a lognormal function, a t low extinctions, with a 
power-law tail, at larger A (e.g., iKainulainen et aL I l2009l . 
see also Ballesteros-Paredes et al. 2011b for the theoretical 
counterpart). The functional form of such iV-pdf is, thus, 



,A A)dA = { Po exp (~ ( ' nA 2°^ m " )2 ^ dA/A if A < A cut , 

Pi(A/A cut )- n dA ifA^A cut) 



where In A max denotes the extinction value at which the log- 
normal is maximum, and In A cut the place where the iV-pdf 
changes from a lognormal to a power-law function. The con- 
tinuity condition at A cu t is: 



Pi = 



Po 
A cu i 



exp 



(In A cut — In A n 
2^2 



(25) 



As in the case of two power-laws ( H3.2J1 . at large extinctions 
(A > A2), the mass, surface and column density must be 
given by the solution of a single powerlaw, eqs. 0, © and 
For smaller extinctions (A < A cu t), however, the solu- 
tions for M and S are: 



/ In (Aut I A max I 

1 v 



S(Ao) = Stot^ Po\/ —o 



V V2a 



(26) 



1 - n 
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Figure 3. Upper panel: 7V-pdf of a hypothetical MC exhibiting 
three peaks. Lower panel: its corresponding £ vs Ao diagram. 
The dotted line denotes the line with slope /9Z , where /3 = 1.67 X 
10 22 cm — 2 mag -1 , the opacity in the K band. 



and 



Figure 4. Upper panel: N-pdl defined by a lognormal function 
at low extinctions and a power-law tail at large extinctions, for 
different values of A cut : 0.2, 0.3, and 0.4. £ = 0.4 and n = 4 
(m = 3) remain fixed. Lower panel: The corresponding £ — Ao 
diagram. Notice that in the power-law regime (A > A cu t), the 
column density is the same for the three cases, since they have 
the same value of n and A cu t, (see eq. J9]l). 



M(A ) = StotfJ,m H p i -paxl^oe" /2 A n 



erf 



erf 



o- 2 - In A cut + In A n 
a 2 — In Ao + In A ma 



V2a 



Pl^cut 

2 - n 



(27) 



while the column density is given by the ratio of eqs. (|27p 
over (|26[) . We now plot the iV-pdf as a function of the extinc- 
tion A (upper panels), and its corresponding mass column 
density E as a function of the extinction threshold Ao (lower 
panels). In Fig. [4] we vary the value of A cu t, the extinction 
value where we the power-law cuts the lognormal. In fig. [5] 
we vary the value of the power-law index n, while in Fig. [6] 
we vary the value of the standard deviation a of the lognor- 
mal function. 

As we can see, again, the column density is basically 
the same in all cases, regardless of the parameter we are 
varying, and large differences occur only when we strongly 
vary the power-law index n. 



3.6 Conclusions from the E — Ao diagrams 

In the previous sections we have constructed different func- 
tional forms of the iV-pdf. Their main features are: (i) they 
peak at the same value (Ai ~ 0.1), (ii) and they decrease 
strongly for larger extinction values. In all cases we have re- 
covered the main two features observed for the column den- 
sity of MCs: at low extinction thresholds, the column density 
is nearly constant, while at large extinction thresholds the 
column density increases linearly with the threshold. Addi- 
tionally, in all cases the values obtained for E at any Ao are 
quite similar, despite their strong differences in functional 
form of the iV-pdf. 

The reason for which very different iV-pdfs produce ap- 
proximately the same values of the observed column density 
at a given threshold, appearing to reproduce the M oc Ft 2 
relation, is a consequence of the definition of E itself, along 
with the fact that the iV-pdfs of MCs decay fast: the col- 
umn density given by the ratio of eq. @ over eq. is the 
average of the extinction coefficient (which is proportional 
to the mean surface density) above the extinction threshold 
Ao. Lets assume now that the iV-pdf has an arbitrary shape 
below its maximum at Ai, and decreases strongly beyond 
that value. For Ao < Ai, the mean column density does 
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Figure 5. Upper panel: Af-pdf defined by a lognormal function 
at low extinctions and a power-law tail at large extinctions, for 
different values of n: 3, 4, 5, and 6. £ = 0.4 and A2 = 0.2 remain 
fixed. Note that the slope that one measures in this plot is m = 
n — 1. Lower panel: The corresponding S — Ao diagram. 

not change substantially with the threshold Ao because the 
mass at a given extinction A, which is proportional to the 
product of the number of points contributing to the JV-pdf 
times A, (dm oc ApdA), will be dominated by the mass at 
the peak A\ . Note that this is valid even if the ./V-pdf were 
flat for Ao < A\. This is the case where the surface density 
is constant with Ao in the £ — Ao diagram. 

As we approach the peak of the iV-pdf, small changes 
in the threshold Ao start removing a substantial number of 
low extinction data points from the iV-pdf. The calculation 
of our mean extinction starts to be dominated by larger val- 
ues of A, making the mean surface density (£} to increase 
rapidly. Once we have passed the maximum (Ao > A max ), 
and the distribution decreases rapidly, the mass will be dom- 
inated by values close -but above- to Ao. This is the range 
in which the (£) varies linearly with the threshold, i.e., 
(£) oc A . 

A similar analysis can be performed if the iV-pdf ex- 
hibits several bumps, as in the case of our three-peaked N- 
pdf (see Fig[3|: below each peak, the E — Ao curve is more 
or less constant because the mass is dominated by the next 
peak. After each peak, the column density grows fast be- 
cause we are removing a substantial number of points at 
low column densities. Once we reach the final peak, the col- 
umn density should increase proportional to Ao because the 
mass is dominated by values slightly larger than Aq. 



Figure 6. Upper panel: N-pdi defined by a lognormal function 
at low extinctions and a power-law tail at large extinctions, for 
different values of a: 0.3, 0.4, 0.5, 0.6, and the unrealistic cases 
with o~ =1 and 1.5. A2 = 0.3 and n = 4 (m = 3) remain fixed. 
Lower panel: The corresponding £ — Ao diagram. As in the case 
of Fig. [4] the column density in the power-law regime (A > A cn t) 
is the same for the three cases, since they have the same value of 
n and .4 cut . 

We emphasi ze that the arguments given bv lKeeell (|l989l l 
and IScald (|l990h . that clouds are not seen at large column 
densities because of optically thick effects, do not apply be- 
cause of the larger dynamic range available in more recent 
observations. Clouds approximately fall along the constant 
column density line just because we are measuring the mean 
value of a property which probability distribution function 
has a peak at nearly the same value for all objects, and falls 
fast enough. 

On the other hand, from the present analysis it becomes 
clear how a molecular cloud must behave in order to depart 
from the Larson's third relation, i.e., in order to have, at 
a given threshold, a substantially different column density 
compared to other clouds: (a) it must have a wider iV-pdf, 
or (b) it must peak at a substantially different extinction 
value A max . The first case can be achieved by having either 
a iV-pdf with a much flatter slope at large column densities, 
or a wider lognormal. Nevertheless, all the iV-pdfs reported 
in the literature fall a factor of ~ 3 orders of magnitude 
when increasing one ord er of magnitude in column density 
l|Kainulainen et ai1l2009l ). making this possibility inapplica- 
ble from the observational point of view. Moreover, it has 
been shown theoretically that, for a supersonic turbulent 
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Figure 7. Schematic T, — Ao diagram for local clouds (shadow re- 
gion), and for some hypothetical clouds (white region). Note that, 
the wider the Af-pdf, the wider the gray region. If the dust opac- 
ity of such hypothetical clouds is the same than that of the local 
clouds, their locus in this diagram should be shifted towards the 
upper right region, as indicated. Here we assu med that A max ~ 1 
for th ese clouds. Observed N-pdi of IRDCs by[Kainulaiiic n"et al.l 
Hon]) discard the possibility of these clouds. 



field with Mach number M, the volume density PDF has a 
standard deviation a p given by 



a logp = ^ln(l + Af 2 /4) 



(28) 



ijQstriker et al.1 1 19991 : IPadoan fc Nordlundl I2002T ). Clearly, 
the dependence on M is weak. If one takes one of the most 
extreme cases known, t he GO. 253+0. 016 mol ecular cloud 
near the galactic center (|Longmore et alj|2012l ). with a ve- 
locity dispersion of 16 km sec -1 , the width of the volume pdf 
will differ by only a factor of two from more typical clouds, 
like Taurus with a dispersion of ~ 2 km sec -1 . But the col- 
umn density PDF by definition must have a smaller width, 
since the column density is the integral of the volume den- 
sity, and thus, density fluctuations are smeared out. Thus, 
one should not expect a wide lognormal iV-pdf even for large 
Mach numbers. 

In the second case, i.e., for hypothetical clouds peaking 
at substantially different extinction values A m&K , if the dust 
opacity /3 is the same for such a cloud, its E — Ao curve will 
be shifted along the line with slope unity, increasing also its 
column density for Ao < ^4 m ax, as shown in Fig. [7) From this 
figure, one may think that for extinction thresholds Ao < 
1, such a cloud will exhibit a much larger column density 
at a given threshold, compared to the Solar Neighborhood 
clouds. Although one may think that this could be the case 
of the infrared dark clouds (IRDCs), which achieve column 
densities substantially larger than local clouds with typical 
values of 10 22 cm" 2 and maximum values up to 10 24 cm -2 



(e.g., iRagan et alj 2011 ). the cu mulative mass functions of 
IRDCs Kainulainen et akl l|201ll . see their Fig. 7) show that 
the mass keeps increasing moving to lower extinction Av ~ 
2, indicating that the peak of their iV-pdf must be at A v ~ 
1 — 2, just as local clouds. 



10 5 

10 4 

1000 
10 5 

g 10 4 

o 

1000 
10 5 

10* 

1000 



ppp 



H — h- 1- 



-I— I- 



CS(l-O) 



i mi} 



H — h- 1- 



— - 



I3 CO(1-0) 




0.01 

size [pc] 



0.1 



Figure 8. Density-size relation for clumps in numerical sim- 
ulations. Upper panel: real space, or position-position-position 
(PPP) space. The clumps are found by thresholding the vol- 
ume density field. Middle and lower panels: observational 
space, position-position-velocity (PPV) space. In this case, 
the clumps are found by thresholding the intensity, which in 
turn is proportional to the column d ensity. Figure taken from 
IBallesteros-Paredes fc Mac Low! [120021) . 



4 THE EFFECT OF THRESHOLDING 
DENSITY 

Imagine now that we were able to measure the volume den- 
sity fields of MCs, and then threshold them in order to mea- 
sure the mean volume density of the clouds. In this case, 
what we would obtain is a nearly constant volume density, 
implying that M oc R 3 . To show that this is the case, in 
Fig. [8] we show the volume density (y axis) against size (x 
axis) for clumps in a numerical simulation of isothermal 
molecular clouds with forced turbule nce, as published by 
IBallesteros-Paredes fc Mac Lowl |2002l . Fig. 9). The details 
of the particular simulation used and on the assumptions 
used for the radiative transfer performed to mimic the line 
emission can be found in that paper. Here we stress that 
the same result can be found in simulations with or without 
self-gravity, with or without magnetic fields, or whether the 
turbulence is forced or decaying. The data points in the up- 
per panel are obtained by thresholding the volume density, 
while the middle and lower panels are obtained by thresh- 
olding the simulated CS(l-O) and the 13 CO (1-0) intensity, 
respectively, which is proportional to the column density. 
The dotted line in these panels has a slope of —1, i.e., its a 
line with constant column density. 

It can be seen from the upper panel that the vol- 
ume density seems to be constant for all cores, when we 
threshold the volume density. When we analyze the obser- 
vational space (middle and lower panels), the points follow 
the straight line with slope —1, implying that the column 
density is constant. In this case, the cores were found by 
thresholding the intensity, which is proportional to the col- 
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umn density. Thus, while we will deduce a relation M oc Ft 3, 
from the upper panel, from the lower two we will deduce 
that the relation is M oc R 2 . 



5 THE INTRA-CLOUD MASS-SIZE RELATION 

We now turn to the intra-cloud mass-size relation, i.e., the 
M — R relation for a single cloud. In this case, Fig. 2 by 
lLombardi etafl (|2010T) show that none of the M — R dia- 
grams exhibit a power-law of the forrr0: 



M oc R a , 



(29) 



Instead, the M — R diagram of all clouds is a curved line with 
the following properties: (a) at the smaller radii, all curves 
exhibit a power-law with slope of 2, and (b) for the larger 
radii, all the curves become flatter, with slopes of the order of 
unity. A natural consequence of these properties is that there 
is a region at intermediate radii (between 0.3 and 3 pc) where 
the curve can be fitted by a power-law with an exponent 
of a ~1.6-1.7. This is the observational result qu oted by 
iKauffmann et al l l|201Ch and lLombardi et al l (|2010l ). 

In order to understand the origin of this behavior, we 
first note that, by construction, the maximum possible value 
of the slope a of the M — R curve is 2: at a given threshold, 
each cloud has a particular size and mass, i.e., each cloud has 
a particular point in the M — R diagram, but all of them fall 
along the same straight line with slope of 2, as a consequence 
of the third Larson's relation. For every cloud, smaller sizes 
and masses are obtained when we increase the extinction 
threshold. To visualize this, pay attention to any pair of 
points with the same symbol in Fig. 1 by lLombardi et al.l 
(2010). In this case, we are moving from one point along a 
line with slope of 2, to another point, with smaller mass and 
size, along another line with slope of 2, but with a larger 
intercept. Then, for a single cloud, their connecting points 
at different thresholds are shifted along lines with slopes 
smaller than two. The limit case occurs when the maximum 
column density of the cloud is reached, and it is in this case 
when q — > 2. 

We now calculate the intra-cloud M — R relation, for 
clouds with different functional forms of their JV-pdfs: log- 
normal, power-law and the combination of these two. Adopt- 
ing the definition 7? = yj S/ir, as in most observational works 
on MCs, we can calculate the mass-size relation for a single 
cloud from the iV-pdf: the mass will be given by eq. (|4]), 
while the size will be proportional to the square root of eq. 
©■ 

Fig. [9] shows the mass-size diagram for four lognormal 
functions: two with sizes similar to those of low mass clouds 
(R = 10 pc, e.g., the Pipe nebula) and two with sizes similar 
to those of large-mass clouds (R = 50 pc, e.g., Orion). Each 
one of these cases have two different values of a: 0.29 and 
0.59, which are the extreme values of t he standard devia- 
tion a of the lognormal fits obtained by iKainulainen et all 
(2009). The solid straight line running across the diagram 



3 We use the letter a to denote the exponent of the hypothetical 
M — R power-law, and a to denote the slope of the curve of the 
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Figure 9. M — R diagram for four lognormal ./V-pdfs. In all cases, 
Amaxi the peak of the distribution, occurs at A max = 0.1. The 
size of the cloud, and the standard deviation of the Gaussian 
are indicated. Lower panel: Same diagram for different power law 
iV-pdfe. 
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Figure 10. M — R. diagram for four power-law ./V-pdfs. In all 
cases, Amax, the peak of the distribution, occurs at Amax 

= 0.1. 

The size of the cloud, and the index of the power-law at large 
extinctions A are as indicated. In all cases, q = 0. Note that the 
index q for low extinctions at large radius. Thus we do not make 
plots for other values of q. The slope of these lines follows eq. I|30|l . 



has a slope of 2. As we can notice, the slope of the mass- 
size relation at s mall radii approaches to 2, as in Fig. 2 by 
lLombardi etal] (|201Ch . 

In Fig. QJJwe show the M — R diagram for four clouds 
with two power-law iV-pdfs. Again, there are two cases with 
size of 10 pc, and other two with 50 pc. The indexes in the 
high extinction wing of the 7V-pdf are n = 3.5 and 6 (m — 
2.5 and 5). In all cases, the index of the low extinction wing 
of the iV-pdf is q = 0. We note that a very steep power-law 
iV-pdf is needed, with an index of n = 6 ( m = 5), in order 
to rep roduce the slope of 1.6 reported by lLombardi et al.l 
(2010) for the Al — R diagram at intermediate radii. 

In Fig. 1111 the mass-size diagram for clouds with a log- 
normal function at low column densities, and a power-law 
regime at large extinctions iV-pdf. The fiducial parameters 
are: A cut — 0.2, n = 4, and a — 0.4. In the upper panel 
we vary the value where the transition between the lognor- 
mal and the power-law occurs: A cut = 0.2, 0.3 and 0.4. In 
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Figure 11. M — R diagram for lognormal+power-law ./V-pdfs. In 
the upper panel, we vary the extinction value where the transition 
between the lognormal and power-law regime occurs, A cu t. In 
the middle panel we vary the slope of the power-law regime, n. 
Finally, in the lower panel, we vary the width of the lognormal 
function, cr. 



the middle panel we vary the slope of the power-law regime: 
n = 6, 5 and 4. Finally, in the lower panel we vary the width 
of the lognormal: a =0.3, 0.4, 0.5, 0.6, 1 and 1.5. 

What these figures show us is that the width of the 
lognormal function shifts the relation: the wider the 7V-pdf, 
the larger the mass at a given radius. Similarly, the slope 
n of the power-law in the iV-pdf determines the slope a of 
the M — R relation: Assuming a po wer- la w iV-pdf, from eqs. 
@ and (JSJl, we can solve for 7? = */ S/tt, and show that the 



Figure 12. n — a diagram, according to eq. I|30|l . where a stands 
for the exponent of the M — R relation, while n stands for the 
exponent of the power-law of the iV-pdf at large extinctions. We 
recall that the slope measured from the iV-pdf when plotted in 
log-log scale is given by m = n — 1, as given in the upper x-axis. 
Notice that value a = 2 is forbidden, according to eq. 1301 



index a in eq. Q29I) will be related to the index n of the high 
extinction wing of the iV-pdf by 

4- 2n 



(30) 



Fig. 1121 shows the n — a diagram, as given by eq. (|30[) . 
We notice that a moves softly from ~ 1.2 to 1.6 for realistic 
values of the power-law (n = 3.5 — 6, or m = 2.5 — 5). 
Note furthermore that value a = 2 is forbidden, according 
to eq. PJ)) . 

All together, our results and the observational results 
present the following puzzle: while most of the observed 
clouds exhibit a power-law distribution iV-pdf for large ex- 
tinctions, why such clouds exhibit an M oc R 2 relation at 
small radii, if the value a — 2 is forbidden for power-laws 
at large extinctions? In principle, this value will be allowed 
only for the few clouds that exhibit lognormal iV-pdfs, as 
Coalsack, California, the Pipe, Lupus V, LDN 1719. 

In order to explain this apparent contradiction, lets as- 
sume that we have a iV-pdf with a shape that changes con- 
tinuously, such that it can be approximated as a sequence 
of power-laws in small intervals of A: 



P (A)dA = pi 



dA for Ae(Ai,A i+1 ), Ai < Ai+i (31) 



We now assume that the iV-pdf decreases rapidly, such 
that the mass at a given threshold is dominated by the low- 
density region (this assumption is fulfilled by the observed 
clouds, since this is one of the conditions needed in order 
to obtain the third Larson relation, as we showed in the 
previous sections). In this case, we can integrate out eqs. §3§ 
and Q only for the small interval in which the power-law 
applies, i.e., 



S(A ) = Stat / p(A)dA 



+ i 



Stot / Pi (A/Ai)- n *dA 



(32) 
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and 



M(A ) 



St t/J,m H P A p(A) dA 



S to tHm H P / A Pi (A/Air ni dA. (33) 



Solving these equations we obtain 

1-71; 



S{Ao) = 



and 



M(A ) 



Stot Aj pi 
1 — m 



Ai+i 
Ai 



Stot ym H P Pi A 2 t 
2-m 



Ai 



(34) 



(35) 



Since we have defined S — tyR 2 , we can write the ratio 
Ai+i/Ai in terms of R from eq. (f34|) : 

(36) 



A i+ i 



1 — Tli „a 
1 + 7ri? 



A» 1 Stot Ai p 
and insert this value in eq. (|35[) to obtain 



M(A ) 



with 



Stot jim H PPiAj 
2 — Hi 



S 'tot AiPi 



7rir+i 



1 (37) 



where 5 < 1. 

Lets now analyze the limit cases. For large radii (low 
extinctions), tyR 2 ~ Stot- Note that the ratio (l — nj)/AiPi ^ 
25: observations by iKainulainen et al.l |2009h show that the 
typical values of the parameters are: 3.5 ^ m 6, 0.1 ^ 
Ai ^ 2, (recall that Av ~ lO^lfe, and our are assumed to 
be in the K band) and Pi ~ 1 for Aj ~ 0.1, and p » 5j 10~ 3 for 
A, ~ 2 (see iV-pdfs in IKainulainen et al.l [20091 ) . Thus, the 
term in parenthesis becomes much larger than unity, and 
then, 



M{A ) oc R 2 



(39) 



for large radii, recovering the solution for the single power- 
law case, eq. (|30[1 . Note that 28 < 2, since 8 < 1. As an 
example, lets take for instance the jV-pdf of Lupus I cloud, 
shown in Fig. 4 by IKainulainen et al.1 |2009). The slope for 
Av ^ 1 is approximately m = 3.3. Thus, n = 4.3 and 28 — 
1.13. This is t he value of the sl ope t hat can be measured 
from Fig. 1 in ILombardi et al.l (|2010ft for the black circles 
that denote Lupus 1. 

Lets now analyze the opposite case. For small radii, 



-tyR 2 « 1 



s . . (40) 

otot Ai pi 

and thus the expansion of eq. (|37p in Taylor series gives the 
relation 



M(A ) = fj.mii P At tyR 



(41) 



i.e., for small radii all the slope of the M — R curves should 
approach to two, the limit case. For intermediate radii, how- 
ever, the solution must be given by eqs. (|37p and fl38|). 

In summary, as we can see from eqs. (|37p ~ (14ip . the intr- 
acloud mass-size relation is completely defined by the shape 



of the iV-pdf, and transitions from a slope of 2, at small 
radii, towards flatter slopes, at large radii. 



6 DISCUSSION 

6.1 The inter-cloud relation 

Early discussions of the importance of the apparent con - 
stancy of the column density of MCs found by |Larso"nl (Il98ll ) 
emphasized the limited colu mn densities probed by then- 
available 12 CO observations (|Kegellll989l ). Since that time, 
other tracers of column density such as near-infrared ex- 
tinction with larger dynamic ra nge, have been used t o de- 
rive similar resul t s. In p articular, ILombardi et al.l (|2010ft and 
IBeaumont et al.1 l|2012ft have found M oc area at differing 
thresholds, though with differing normalization. 

In the first sections of this work we showed how a wide 
variety of iV-pdfs can reproduce the inter-cloud Larson's 3rd 
law, i.e., can produce M oc area relations for multiple clouds. 
Furthermore, we showed that if one could in principle thresh- 
old the volume density, one would get a M oc volume rela- 
tion. This strongly suggests that these types of relations, 
rather than universal laws, are mostly an artifact of the 
thresholding process. 

Operationally, cloud masses are obtained by integrating 
or summing column densities over the cloud area. It is thus 
not surprising that cloud masses are observed to correlate 
with areas, specially, if most of the mass is at low densities, 
i.e., occupying most of the area. If then the JV-pdf is rapidly 
decreasing with increasing column density, this means that 
most of the column density and thus most of the mass is near 
the threshold. Thus one automatically gets a cloud mass 
linearly proportional to its area. Similarly, if one were able 
to measure densities directly, the same principle results in 
mass linearly proportional to volume. 

The significance of the mass-area relationship then 
reduces to understanding why the mean or peak of the 
iV-pdf seems not to vary much from region to region 
ijBeaumont et all |2012| ). and why the iV-pdf has a rapid 
decrease with increasing column density in the thresholds 
of interest. We argue that there are good physical reasons 
for this. The first is that, at least in the Solar Neighbor- 
hood, an extinction of Av ~ 1 — 2, equivalent to a hydrogen 
column density of N ~10 21 cm -2 , is approximately what 
is needed to shield CO from the dissociating photons of 
the interstellar radiation field (Ivan Dishoeck fc Blacklll988l : 
Ivan Dishoeck fc Blakdl 1998ft . It is difficult to set thresholds 
for molecular clouds below this average value (which corre- 
sponds to A ~ 0.1 at K-band, the lowest level considered by 
ILombardi et al.ll2010ft . since as it would not be a CO cloud 
and these studies are focused towards MCs. 

Secondly, we argue that most of the mass of clouds is at 
low density, because the shielding value is roughly the col- 
umn den sity at which one expects self-gravity to becom e im- 
portant (|Franco fc Coxlll986l: iHartmann et alj|20qift. Nu - 
merical simulations (e.g.. iBallesteros-Paredes et alj|2011bft 
indicate that the slope of the iV-pdf at large column densi- 
ties iV is steeper than 3, as a result of gravitational collapse. 
More specifically, in non-spherical, non-uniform geometries 
gravitational collapse can be highly non-li n ear as a function 
of position (jBurkert fc Hartmannl 120041 : IHartmann et al.1 
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2012), such that the dense material is a small fraction of 

the total cloud mass. 

iBeaumont et al.l l|2012h concluded that the "univer- 
sality" of Larson's third relation is a result of having 
the dispersion in the mean value of the column den- 
sity be smaller than the dispersion in cloud areas. Our 
argument is that this occurs because the large column 
density regions of clouds have a small filling factor as 
a consequence of the gravitational contraction and col- 
lapse, producing the m ass oc area result. We agree with 
IBeaumont et all |2012l 'l that therefore the iV-pdf itself is 
much more instructive about cloud structure than the mass- 
area "law", and that the power-law or other high-column 
density "tails" on iV-pdfs i s a resu lt of gravitational collapse 
( Balle steros-Paredes et ail I2011b1 ). This is consistent with 
our previous ar guments that mod ern versions of Larson's 
second relation jHever et al.ll2009f) imply global (and l ocal) 
gravitational collapse 1 Ballesteros-Paredes et al.ll2011aT ). not 
Virial equilibrium. 



6.2 The intra-cloud relation 



Although the original mass-size relation by iLarsonl (|l98ll ) 
M oc R 2 was obtained for a set of clouds, i.e., in terms 
of the inter-cloud relation, different authors had noticed 
that this relation does not holds for a single cloud. The 
typical values of t he exponent used to be smaller, typi- 
cally M oc r 1 ' 6-1 - 7 (|Kauffmann et alj|201(il : lLombardi et all 
2010), altho ugh it should be recognized that the diagrams 
presented by lLombardi et all (|2010T l show that the relation 
is a curved line, rather than a power-law. The results pre- 
sented in the present paper shows that the iV-pdf of the 
cloud determines the slope of the M — R relation, and it 
goes from 2 as a limit case for small radii, and flattens for 
larger radii. 



7 CONCLUSIONS 

By using simple functional forms of the column density prob- 
ability distribution function (iV-pdf), we have explored the 
significance of the third Larson scaling relation for molecu- 
lar clouds. We show that for a set of clouds, this relation is 
an artifact of thresholding the surface density, as long as the 
form of iV-pdf peaks near or below the threshold used, and 
that the mass is a steeply decreasing function of increasing 
column density. We argue that the physical reasons for these 
features of the column density distribution are that there is 
a minimum column density to shield CO, and that this min- 
imum value is not far from that at which the molecular gas 
becomes Jeans unstable, resulting in cloud collapse to form 
stars. 

We also showed that for a single cloud, the slope of 
the mass-size relation approaches to 2 at small radii, and 
flattens for larger radii, in agreement with observations. The 
detailed shape of this curve depends on the detailed shape 
of the iV-pdf. 
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